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The performance of Navier-Stokes codes are influenced by several phenomena. For 
example, the robustness of the code may be compromised by the lack of grid resolution, 
by a need for more precise initial conditions or because all or part of the flowfield 
lies outside the flow regime in which the algorithm converges efficiently. A primary 
example of the latter effect is the presence of extended low Mach number and/or low 
Reynolds number regions which cause convergence deterioration of time marching 
algorithms. Recent research into this problem by several workers including the present 
authors has largely negated this difficulty through the introduction of time-derivative 
preconditioning. In the present paper, we employ the preconditioned algorithm to 
address convergence difficulties arising from sensitivity to grid stretching and high aspect 
ratio grid cells. 

Strong grid stretching is particularly characteristic of turbulent flow calculations 
where the grid must be refined very tightly in the dimension normal to the wall, without 
a similar refinement in the tangential direction. High aspect ratio grid cells also arise in 
problems that involve high aspect ratio domains such as combustor coolant channels. 

In both situations, the high aspect ratio cells can lead to extreme deterioration in 
convergence. It is the purpose of the present paper to address the reasons for this 
adverse response to grid stretching and to suggest methods for enhancing convergence 
under such circumstances. 

Numerical algorithms typically possess a maximum allowable or optimum value 
for the time step size, expressed in non-dimensional terms as a CFL number or von- 
Neumann number (VNN). In the presence of high aspect ratio cells, the smallest 
dimension of the grid cell controls the time step size causing it to be extremely small, 
which in turn results in the deterioration of convergence behaviour. For explicit schemes, 
this time step limitation cannot be exceeded without violating stability restrictions of 
the scheme. Oh the other hand, for implicit schemes, which are typically unconditionally 
stable, there appears to be room for improvement through careful tailoring of the time- 
step definition based on results of linear stability analyses. In the present paper, we 
focus on the central-differenced alternating direction implicit (ADI) scheme. The 
understanding garnered from this analyses can then be applied to other implicit 
schemes. 

In order to systematically study the effects of aspect ratio and the methods of 
mitigating the associated problems, we use a two pronged approach. We use stability 
analyses as a tool for predicting numerical convergence behavior and numerical 
experiments on simple model problems to verify predicted trends. Based on these 
analyses, we determine that efficient convergence may be obtained at all aspect ratios 
by getting a combination of things right. Primary among these are the proper definition 
of the time step size, proper selection of viscous preconditioner and the precise treatment 
of boundary conditions. These algorithmic improvements are then applied to a variety of 
test cases to demonstrate uniform convergence at all aspect ratios. 
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Philosophy of Grid Aspect Ratio Study 


• Assessment of High Aspect Ratio Problem 

— Disparate propagation speeds in X and Y 

• Stability Theory 

— Scalar Convection-Diffusion Equation 

— Euler Equations 

— Navier- Stokes Equations 

• Numerical Convergence Studies 

— Simple Model Problems 

— Realistic Flow Problems 

• Improved Algorithm to Provide Aspect Ratio Control 

— Precise Time-Step Definition 

— Viscous Preconditioning 

— Boundary Condition Implementation 
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The Navier-Stokes Equations 


r dQv dE dF 
dt dx dy 


H + L(Q V ) 


• Solution Vector 


Qv = ( P, u, v, T ) T 

• Preconditioning Matrix 
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• Parameter e 

— Activates Inviscid and Viscous Preconditioning 

— Value Depends on Local Mach Number and Cell 
Reynolds Number 
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Numerical Solution Procedure 


Central-Differenced ADI Algorithm 
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• Approximate factorization errors control convergence 
behavior. 

• Optimum CFL u+c is typically between 1 and 10. 

— Other inviscid and viscous time scales are opti- 
mized by the preconditioning matrix. 
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Time-Step Definition 


• Local Time-Stepping or Constant CFL Condition 

Max ( CFL X , CFL y ) = CFL 


• For high aspect ratios, CFL X and CFL y become dis- 
parate 


CFLy = CFL , CFL X = CFL/AR 
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Euler Stability Analysis 


• Aspect Ratio ( AR ) of 100 
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New Time-Step Definition 


• Conclusions from Stability Analysis: 

— Min-CFL Preferable to Ma x-CFL 
— Efficient Convergence at all AR 

• Minimum-CFL Definition 

Min ( CFL X , CFL y ) = CFL 

• For high aspect ratios, 

CFLy = CFL * AR , CFL X = CFL 
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Implementation of Boundary Conditions 


• Extrapolation vs Characteristic 

— Both work well for small CFL ' s 

— Characteristics usually superior at high CFL' s 

• Proper MOC Implementation: 

— Implicit procedures 

— Boundary conditions applied before approximate fac- 
torization 

— Consistent order of accuracy: LHS / RHS 
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No. of Iterations to Machine Accuracy 


High Aspect Ratio Convergence — Inviscid Duct 
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Navier-Stokes Analysis 


• Parameter e controls low Re number convergence 

— Viscous terms limit time step at high AR 

— e chosen to optimize inviscid and viscous modes simul- 
taneously 

• Obvious choice: e = f (Ma x-CFL, Ma x-VNN) 

• Scalar Stability Results: e = / (Min-CFL, Min -VNN) 
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Navier-Stokes Analysis 


• Min-CFL, Min-V'iVjV Stability Result 
AR = 1000 



VNN X = 1, VNN y = 1 x 10 6 
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Navier- Stokes Analysis 


Conclusions from Stability Results: 

— Vector system different from scalar equation 

— Approximate factorization error CFL X * VNN y limits 
convergence 

Viscous preconditioner, e = / (Min-CFL,M&x-VNN) 

— Maintains Min-C'FL for ‘inviscid’ modes 

— Uses traditional Max-WiV definition for ‘viscous’ modes 



Navier- Stokes Analysis 


• Min-CFL, Mbx-VNN Stability Result 
AR = 1000 



to x 


CFL X = 1, CFL y = 1000 
VNN X = 1 x lO" 6 , VNN y = 1 
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No. of Iterations to Machine Accuracy 


High Aspect Ratio Convergence — Viscous Duct 



Grid Aspect Ratio 
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RADIAL COORDINATE 
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Hydrogen/ Oxygen Shear-Layer 
Convergence 



No. of Iterations 
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L 2 Norm of the Residual 


High Reynolds Number Boundary Layer 
Convergence — Stretched Grid 



No. of Iterations 
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High Reynolds Number Boundary Layer 

Blasius Solution 



u/U 
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High Reynolds Number Boundary Layer 
Turbulent Velocity Profile 
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Turbulent Nozzle Computation 
Stretched Grid and Solution 
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Ij 2 Norm of the Residual 


Turbulent Nozzle Computation 
Convergence 



No. of Iterations 
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Heat Flux, W/m 


Turbulent Nozzle Computation 
Wall Heat Flux 
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Conclusions — High Aspect Ratio Study 


• High Aspect Ratio Analysis 

— Stability Theory 

— Numerical Convergence Studies 

• Convergence Control: 

— Min -CFL Time Step 

— Max-WiV Viscous Preconditioner 

— Correct implementation of boundary conditions 

• Uniform Convergence Demonstrated for All AR’s 

— Above issues addressed in combination 

— Efficient convergence for variety of test cases 
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Conclusions (Contd.) 


• Present results are for two-dimensional central-differenced 
ADI scheme 

• Explicit Schemes: 

— Optimum time step causes poor convergence at high 
A, R’s 

• Upwind Schemes Also Suffer at High AH’s 

— Present improvements may be incorporated 

— Rich variety of approximate factorization methods 

• Three-dimensional computations: 

— ADI scheme is conditionally stable 

— Two kinds of high aspect ratio grids 

— Algorithmic improvements appear promising 
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